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Abstract 

We discuss production of nonphotonic electrons in proton-proton scattering at RHIC. The dis- 
tributions in rapidity and transverse momentum of charm and bottom quarks / antiquarks are cal- 
culated in the fct-factorization approach. We use different unintegrated gluon distributions from 
the literature. The hadronization of heavy quarks is done by means of Peterson and Braaten et al. 
fragmentation functions. The semileptonic decay functions are found by fitting recent semileptonic 
data obtained by the CLEO and BABAR collaborations. We get good description of the data at 
large transverse momenta of electrons and find a missing strength concentrated at small transverse 
momenta of electrons. Plausible missing mechanisms are discussed. 
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I. INTRODUCTION 



Recently the PHENIX and STAR collaborations has measured transverse momentum 
distribution of so-called nonphotonic electrons [H, Hj. It is believed that the dominant contri- 
bution to the nonphotonic electrons/positrons comes from the semileptonic decays of charm 
and beauty mesons. These processes have three subsequent stages. First cc or bb quarks are 
produced. The dominant mechanisms being gluon-gluon fusion and quark-antiquark anni- 
hilation. Next the heavy quarks/antiquarks are turned to heavy charmed mesons D,D* or 
B,B*. The vector D* and B* mesons decay strongly producing D and B mesons. Finally 
the heavy pseudoscalar mesons decay semileptonically producing electrons/positrons. 

The inclusive heavy quark/antiquark production can presently be calculated at Fixed- 
Order plus Next-to-Leading-Log (FONLL) level j^. The predictions for electron spectra in 
proton-proton collisions at RHIC can be found in Ref . N] . An alternative approach for inclu- 
sive heavy quark production is /c^-factorization [1, 0, 0, S, [s, 13, 111 • In this approach emission 
of gluons (see FigH]) is encoded in so-called unintegrated gluon distributions (UGDFs). The 
latter approach is very efficient for description of QQ correlations (l2| . 

The hadronization of heavy quarks is usually done with the help of fragmentation func- 
tions. The Peterson fragmentation functions are often used in this context [3]. The pa- 
rameters of the Peterson fragmentation functions are adjusted to e~^e~ or pp production of 
heavy mesons. Another perturbative fragmentation model has been proposed in Ref. [3] 
(BCFY). 

The last ingredient are semileptonic deca^ of heavy mesons. Until recently this com- 



ponent was treated by modeling the decay [15|, [ly, ll7|. Only recently the CLEO [18j and 
BABAR fl^ collaborations has measured very precisely the spectrum of electrons/positrons 
coming from the decays of D and B mesons, respectively. This is done by producing res- 
onances: \Ef(3770) which decays into D and D mesons (CLEO) and T(4S') which decays 
into B and B mesons (BABAR). In both cases the heavy mesons are almost at rest, so in 
practice one measures the meson rest frame distributions of electrons/positrons. 

In the present analysis we shall apply the fcrfactorization approach. At relatively low 
RHIC energies rather intermediate x- values become relevant. The Kwiecinski unintegrated 
gluon (parton) distributions seem relevant in this case [20\ . We shall use also Ivanov-Nikolaev 
distributions which were fitted to deep-inelastic HERA data including intermediate-x region 
(2l| . We shall use both Peterson and BCFY fragmentation functions. The electron/positron 
decay functions will be fitted to the recent CLEO and BABAR data. 



II. FORMALISM 



Let us consider the reaction hi + h2-^Q + Q + X, where Q and Q are heavy quark and 
heavy antiquark, respectively. In the leading-order (LO) approximation within the collinear 
approach the quadruply differential cross section in the rapidity of Q (yi), in the rapidity of 
Q (7/2) and the transverse momentum of one of them [pt) can be written as 

^^^^ = -mi-^^^') \M-? . (2.1) 

Above, pj(xi,/i^) and Pj(x2,/i^) are the familiar (integrated) parton distributions in hadron 
hi and /i2, respectively. There are two types of the LO 2^2 subprocesses which en- 
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ter Eg. (12. II) : gg QQ and qq QQ. The first mechanism dominates at large ener- 
gies and the second one near the threshold. The parton distributions are evaluated at: 



Xi = ^ (exp(?/i) + exp(?/2)), 2:2 = ^ (exp(-yi) + exp(-?/2)), where rrit = ^pl + The 
formulae for matrix element squared averaged over the initial and summed over the final 



spin polarizations can be found e.g. in Ref . |22 



If one allows for transverse momenta of the initial partons, the sum of transverse momenta 
of the final Q and Q no longer cancels. Formula (12.1 p can be easily generalized if one allows 
for the initial parton transverse momenta. Then 

da sr^ f d^Kitd^K2t 1 



6'^ + K2,t -Pl,t -P2,t) J J^j{x2,nl^t) ) (2-2) 



|2 



dyidy2d'^Pi^td'^P2,t ^ J n n 16tt'^{xiX2s)^ 



where now Ti{x\, /tf J and Tj{x2-, fi^l.t) the so-called unintegrated gluon (parton) distribu- 
tions ^ . The extra integration is over transverse momenta of the initial partons. The two ex- 
tra factors l/vr attached to the integration over rf^Ki^j and (i^K2,t instead over dn^ ^ and dn^ ^ as 
in the conventional relation between the unintegrated (JF) and the integrated ((7) parton dis- 
tributions. The two-dimensional Dirac delta function assures momentum conservation. Now 
the unintegrated parton distributions must be evaluated at: X\ = exp(|/i) + exp(?/2), 

X2 = ■-^exp(— yi) + ^^exp(— 2/2), where rrii^t = y^P^7+^^- ■'■^ general, the matrix ele- 
ment must be calculated for initial off-shell partons. The corresponding formulae for initial 
gluons were calculated in la, \& (see also l?!). It is easy to check [12| that in the limit k1 0, 



/«2 ^ the off-shell matrix elements converge to the on-shell ones. 
Introducing new variables: 



Qt = K,l,t + K,2,t , 

Qt = f^i,t - f^2,t (2.3) 



we can write: 



daii f ,0 1 



dyidy2d'^Pi^td'^P2,t J ^ 4tt2 i67r2(xiX2s)2 

J^i{xi, Tj{x2, K2,t) • (2.4) 

This formula is very useful to study correlations between the produced heavy quark Q and 
heavy antiquark Q (l2| . 
For example 



dp 



^ — = f d(j)id(j)2 Pi,tP2,t [ dyidy2 [ d'^qt \M 

i,tdp2,t J J J l<6ii\xiX2sY 



4- \\ jd'i^- ^1.*^^.* jdy^dy- /^'^* ^^ i&.^{x,x2sf "-^^ 
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J^i{xi, k{,) Tj{x2, KoJ . (2.5) 



^ In this paper we shall use the following convention of unintegrated gluon distributions: ^{x^ K^)dK^ 
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In the last equation we have introduced 0_ = 0i — 02, where 0_ G (-27r, 27r). The factor 
4 vr comes from the integration over 0+ = 0i + 02- The first factor 1/2 comes from the 
jacobian transformation while the second factor 1/2 takes into account an extra extension 
of the domain when using 0_|_ and 0_ instead of 0i and 02- 

At the Tevatron and LHC energies the contribution of the gg QQ subrocess is more 
than an order of magnitude larger than its counterpart for the qq — >■ QQ subprocess. At 
RHIC energy the relative contribution of the quark- ant iquark annihilation is somewhat 
bigger. Therefore in the following we shall take into account not only gluon-gluon fusion 
process i.e. i=0 and j=0 but also the quark-antiquark annihilation mechanism. 

The purely perturbative^ /cf-factorization formalism to hih2 QQ applies if /ti(,K2i > 
Kq. The choice of Kq is to a large extent arbitrary. In Refs. \^ a rather large was chosen and 
the space /^i ^ x ^2 1 "was subdivided into four disjoint regions. For example the contribution 
when both kI ^ and small was replaced by the leading-order collinear cross section. 

Such an approach assures that a^Q > cr^^ (collinear LO) by construction. 
It is rather obvious that the resulting cross section strongly depends on the choice of Kq 
which makes the procedure a bit arbitrary. Our philosophy here is different. Many models of 
UGDF in the literature treat the soft region explicitly. Therefore we use the fcf-factorization 
formula everywhere on the nf^ x plane. 

The production of electrons/positrons is a multi-step process. The whole procedure of 
electron/positron production can be written in the following schematic way: 

" Dq^d ® fn^e , (2.6) 



dyd^p dyd^p 

where the symbol ® denotes a generic convolution. The first term responsible for production 
of heavy quarks/ ant iquarks is calculated in the /c^-factorization approach. Some details 
were already discussed above. Next step is the process of formation of heavy mesons. We 
follow a phenomenological approach and take Peterson and Braaten et al. fragmentation 



functions with parameters from the literature (see e.g. [27!]). The electron decay function 
should account for the proper branching fractions. The latter are known experimentally (see 
e.g. 18, 27|). These functions can in principle be calculated 15, 13]. This introduces. 



however, some model uncertainties and requires inclusion of all final state channels explicitly. 
An alternative is to use experimental input. The decay functions have been measured only 



recently [18|, [IQj. How to use the recent experimental information will be discussed in the 
next section. 



III. RESULTS 

In principle, the semileptonic decays can be modeled (see e.g. 15, [l^, [13]). Since there 
are many decay channels with different number of particles this is not an easy task. In 
our approach we take less ambitious but more pragmatic approach. In Fig]2] we show our 
pur ely mathematical fit to not absolutely normalized data of the CLEO (l8| and BABAR 
[19] collaborations. We find a good fit with: 

fcLEo{p) = 12.55(p + 0.02)2-55(0.98 -p)2-^^ (3.1) 



when both UGDFs are generated perturbatively 
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for the CLEO data [is) and 



fBABARip) = (126.16 + 14293.09 exp(-2.241n(2.51 - 0.97pf) 

(-41.79 + 42.78 exp(-0.5(|p - 1.27|)/1. 8)^-^^) ^ ' ^ 



for the BABAR data 1^. In these purely numerical parametrizations p must be taken in 
GeV. 

After renormahzing to experimental branching fractions for D ^ e (about 10 % ^ 
and 5 ^ e (10.36 ± 0.06(stat.) ± 0.23(syst.) % [19] ) we shall use them to generate 
electrons/positrons in the rest frame of the decaying D and B mesons in a Monte Carlo 
approach. We shall neglect a small effect of the non-zero motion of the D mesons in the 
case of the CLEO experiment and of the B mesons in the case of the BABAR experiment. 
This effect is completely negligible. 

For illustration of the whole procedure in Fig|3] we show as an example two-dimensional 
distributions in rapidity and transverse momentum for charm quarks, D mesons and elec- 
trons from the decay of D mesons. Both fragmentation and semileptonic decays cause 
degradation of transverse momentum. On average pt,e < Pt,D < Pt,c- The spectra of elec- 
trons are much softer than initial spectra of charm quarks. On the other hand the distri- 
butions in rapidity of electrons are much broader than the corresponding distributions of 
quarks/ antiquarks. 

Now we shall concentrate on invariant cross section as a function of electron/positron 
transverse momentum. Such distributions have been measured recently by the STAR and 
PHENIX collaboration at RHIC [H, 3- In Figg] and Fig|5] we show results obtained with 
Kwiecihski UGDF M and different combinations of factorization and renormalization scales 



as well as for different fragmentation functions (Peterson and BCFY). The differences be- 
tween results obtained with different combinations quantify theoretical uncertainties. Simi- 
larly as for the standard collinear approach [J] one gets uncertainties of the order of a factor 
2. We show individual contributions of electrons/positrons initiated by c/c or b/b. The con- 
tribution of the c/c (dashed) dominates at low transverse momenta of electrons/positrons. 
At transverse momenta of the order of 4 - 5 GeV the both contributions become comparable. 
We obtain rough agreement for large transverse momenta. Similarly as for the higher-order 
collinear approach ^ there is a missing strenght at lower transverse momenta. A better 
agreement is obtained with renormalization scale taken as transverse momentum of the ini- 
tial gluon(s). There are two strong coupling constant in the considered order. In practice 
we take as{kl^)as{k2t), i.e. different argument for each running coupling constant. This is 
rather a standard prescription used in fct-factorization approach (see e.g. [sl, [13]) although 
does not have a deep theoretical foundation. In the latter case to avoid Landau pole we use 



analytic prescription of Shirkov and Solovtsov [231 . 

The situation for the Kwiecihski UGDF is summarized in FiglH] where we have shown 
uncertainty band of our theoretical calculation. The upper curves are for = kf and 

= Atjiq and the lower curves are for = ArriQ and fi'j^ = ^toq. Up to now we have pre- 
sented only the PHENIX collaboration data which span a broader range of lepton transverse 



The branching fraction for different species of D mesons is different: 
m.{D+ e+t/eX) = f6.f3±0.20(stat.)±0.33(syst.)%,BR(i:>o ^ e+i/e^)=6.46±0.f 7(stat.)±0.13(syst.)%) 



18j. Because the shapes of positron spectra for both decays are identical within error bars we can take 



the average value and simplify the calculation. 
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momenta. In Figl6]we show also the STAR collaboration data. The experimental results of 
both groups are not completely consistent. In the interval 3 GeV < pt(lepton) < 6 GeV, the 
STAR data points are somewhat higher than the PHENIX data points. This disagreement 
needs further explanation. Our results are roughly consistent with both experimental sets 
at large pt(lepton). There is a missing strenght at small transverse momenta where only the 
PHENIX collaboration data exist. This will be discussed further in the following. 

In FigE] we show results obtained with Ivanov-Nikolaev UGDF. Although there is some 
improvement at low transverse momenta, the cross section for larger transverse momenta 
exceed the experimental data. 

It is not clear for the moment what is the missing strength. Up to now we have included 
only gluon-gluon fusion which is known to be dominant contribution at large center-of-mass 
energies (Tevatron, LHC). At RHIC energies the typical longitudinal momentum fractions 
of gluons are still not too small Xi,X2 ~ 0.01 and the contribution of the quark- ant iquark 
annihilation may be not negligible. Therefore in the following we shall include also quark- 
antiquark annihilation process. Those processes can be included in a similar way in the 
formalism of unintegrated parton distributions. The corresponding diagram is shown in 



FiglHl The Kwiecihski formalism [20] allows to calculate unintegrated quark/antiquark dis- 
tribution in the same framework as unintegrated gluon distributions. In Fig 19] we present the 
contribution of quark-antiquark annihilation qq — >■ cc (dash-dotted line). This contribution 
is similar in size to the gg bb contribution. The contribution of qq bb is negligible and 
is not shown here. 

Study of nonphotonic and hadron correlations allows to "extract" a fractional con- 
tribution of the bottom mesons B/{D + S) as a function of electron/positron transverse 
momentum |2^]. Recently the STAR collaboration has extended the measurment of the rel- 
ative B contribution to electron/positron transverse momenta ~ 10 GeV [2^. In FigJTOl and 
[TT] we present our results for different unintegrated gluon distributions and different frag- 
mentation functions. There is a strong dependence on the factorization and renormalization 
scale in the case of the Kwieciriski unintegrated gluon distributions. A better agreement is 
obtained with the Peterson fragmentation functions. The separation into charm and bot- 
tom contributions is very important in the context of identifying the missing strenght. A 
new correlation method was proposed recently to identify and separate charm and bottom 
production on a statistical basis [25| . The method was tested using known event generators. 
An alternative method of extracting the relative B contribution from azimuthal angular 
correlations of nonphotonic electrons and Dq mesons was proposed [1^. One can hope that 
application of the new methods will help in disantagling the contributions better. 



IV. DISCUSSION OF THE RESULTS 

We have calculated inclusive spectra of nonphotonic electrons/positrons for RHIC energy 
in the framework of the fcrfactorization. We have concentrated on the dominant gluon-gluon 
fusion mechanism and used two recent unintegrated gluon distribution functions from the 
literature. Special emphasis was devoted to the Kwiecihski unintegrated gluon (parton) 
distributions. In this formalism, using unintegrated quark and antiquark distributions, one 
can calculate in addition the quark-antiquark annihilation process including transverse mo- 
menta of initial partons (quarks/antiquarks). In addition, we have used unintegrated gluon 
distributions constructed by Ivanov and Nikolaev to describe deep-inelastic data measured 
at HERA. 
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When calculating spectra of charmed {D, D*) and bottom {B, B*) mesons we have 
used Peterson and Braaten et al. fragmentation functions with model parameters from the 
literature. There are no big differences between results obtained with both fragmentation 
functions. 

A very important ingredient, which influences the final spectra, is the distribution of 
electrons/positrons from the decay of D and B mesons. Here we have used recent results 
of the CLEO and BABAR collaborations. The momentum spectra of electrons/positrons 
from the decays of D and B mesons produced in the e^e~ collisions were used in the present 
calculation to generate distribution of electrons/positrons coming from the decays of D and 
B mesons produced in the hadronic reactions. This way we have avoided all uncertainties 
associated with modeling semileptonic decays of mesons. 

We have compared results obtained in our approach with experimental data measured 
recently by the PHENIX collaboration at RHIC. We get a reasonable description of the 
data at large transverse momenta of electrons/positrons. Similarly as for the higher-order 
collinear approach there is a missing strength at lower transverse momenta. 

Up to now there is no clear explanation of the enhanced production of electrons/positrons 
at low transverse momenta. The uncertainties related to the choice of factorization and 
renormalization scale seems to be insufficient. There can be several reasons of the unex- 
plained strength at low transverse momenta. 

The fcf-factorization approach includes many higher-order contributions which are em- 
bodied in unintegrated gluon (parton) distributions. Some higher-order contributions are 
definitly not included. A simple and transparent example are emissions of gluons of the 
heavy quarks/ ant iquarks. This contribution can be estimated in the standard collinear 
approach. This effect is, however, not limited to low transverse momenta. 

It is commonly assumed that D / D mesons are produced via fragmentation of c/c quarks. 
However, at lower energies (fixed target experiments) an asymmetry between different species 
of D mesons have been observed [28|. This asymmetry can be due to fragmentation of light 
(u,d,s) quarks/antiquarks j29| (g D{qc)c or q D{qc)c) or meson cloud effects sO]- The 



asymmetry increases with rapidity (or Feynman xp). This makes questionable the common 
assumption that D mesons are produced exclusively via fragmentation of c quarks. In this 
context, it would be very useful to analyze electronic spectra at larger rapidities. If these 
mechanisms are responsible for the missing strength then the discrepancy there would be 
even larger. In the moment only muons were measured at forward rapidities [3l[ and there 
seems to be some enhancement, although systematic error bars are rather large. 

The results of the PHENIX collaboration were obtained by subtraction of several com- 
ponents, including decays of vector mesons, so-called Dalitz decays, i^^es decays and other 
mechanisms. All of them "are concentrated" at low transverse momenta [2]. Only a sketch 
of the subtruction procedure was presented j32|. The details of the subtraction are not pre- 
sented m extenso. It is therefore not clear to us how reliable such subtraction is. In addition, 
there are several mechanisms which were not included. These are Drell-Yan processes, pro- 
cesses initiated by two photons (they are expected to be concentrated at low transverse 
momenta) and several other exclusive processes never calculated in the literature. It seems 
therefore difficult to draw definite conclusions before cross section for all these processes 
is evaluated. We leave such calculations for separate detailed studies. In principle, also 



^ There is a substantial fragmentation of light quarks (q 7^ s) in the case of kaon production. Such a 
contribution for D mesons is therefore also not excluded. 
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analysis of coincidence spectra, e.g. in invariant mass of the dilepton pair Mgg, could help 
to pin down the missing mechanisms. 
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FIGURES 




FIG. 1: A basic diagram relevant for gluon-gluon fusion in /ct-factorization. 
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FIG. 7: Transverse momentum distribution of electrons/positrons obtained with Ivanov-Nikolaev 
UGDF and Peterson (left panel) and BCFY (right panel) fragmentation functions. 
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FIG. 8: A diagram for quark-antiquark annihilation. 
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FIG. 9: Transverse momentum distribution of electrons/positrons with Kwiecinski UPDFs. The 
dash-dotted line corresponds to the qq bb contribution. 
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FIG. 10: The fraction of the B decays for the Kwiecihski UGDF and different combinations of the 
factorization and renormahzation scales. 
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FIG. 11: The fraction of the B decays for the Ivanov-Nikolaev UGDF and two different models of 
the fragmentation. In this calculation = 4mQ. 
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